General Disclaimer 


One or more of the Following Statements may affect this Document 


• This document has been reproduced from the best copy furnished by the 
organizational source. It is being released in the interest of making available as 
much information as possible. 


• This document may contain data, which exceeds the sheet parameters. It was 
furnished in this condition by the organizational source and is the best copy 
available. 


• This document may contain tone-on-tone or color graphs, charts and/or pictures, 
which have been reproduced in black and white. 


• This document is paginated as submitted by the original source. 


• Portions of this document are not fully legible due to the historical nature of some 
of the material. However, it is the best reproduction available from the original 
submission. 


Produced by the NASA Center for Aerospace Information (CASI) 





(NASA-CF-1452C1) SSFi STUDIES Of THE SLOrf N77-2J966 

TUflBLIVii 0? VANADYL SFIS P5CBES It. N2HATIC 
LIQUID CRYSTALS Interim Progress Report 

(Texas Univ.) 15 p HC A02/MF 101 CSCL 20L Uncias 

vJ 3/7 b J012t> 


Interim Progress Report 
Grant NSG 1316 

ESR Studies of the Slow Tumbling of 
Vanadyl Spin Probes In Nematic Liquid Crystals 


by 


Dr. Michael P. Eastman 
Dr. Gerald V. Bruno 
Dr. Juan 0. Lawson 


The University of Texas at El Paso 
El Paso, Texas 79968 





% 






MAY 1377 v 

RECEIVED 
nasa sn Fwaurc 
INPUT BRANCH W 


Abstract 


ESR line shapes that are appropriate for slowly tumbling 
vanadyl spin probes in viscous nematic liquid crystals have been 
calculated by the stochastic Liouville method. Because of the 
symmetry possessed by vanadyl, the analysis and interpretation 
of these line shapes is simplified considerably. Spectral line 
shapes agreed well with experimental spectra of VOAcAc in the 
nematic liquid crystal Phase V and BEPC. Deviations from 
Brownian rotational diffusion were noted. A slowly fluctuating 
torque analysis yielded good agreement with the experimental 
spectra. 



Introduction 


The application of paramagnetic spin labels in a variety of 

experiments has resulted in obtaining a wide range of structural 

and dynamical information. ^ In particular, the class of complexes 

2+ 

containing the vanadyl ion VO is being used increasingly as a 
2 

spin probe. In a recent paper, we have applied the stochastic 

Liouville method to calculate ESR line shapes in the slow tumbling 

3 

region for vanadyl complexes in an isotropic media. This analysis 

allows one to extract quantitative information for the rotational 

_ 1 0 _ 7 

correlation times ( x )in the region 10 -10 sec/rad. These slow 

tumbling ESR spectra occur when the vanadyl ion is rigidly attached 

o 

to a macromolecule (radius 5-50 A) or when a vanadyl complex is 
dissolved in a viscous medium. 

Currently, there is much interest in the dynamical properties 

of liquid crystals. The purpose of this paper is to study spin 

2 + 

probes in viscous nematic liquid crystals. Because many VO com- 
plexes have very nearly axially symmetric magnetic parameters which 
is the same degree of symmetry possessed by nematic liquid crystals, 
the subsequent theoretical analysis is simplified and the results 

are easier to interpret than spin probes with less symmetry. Previous 

2 + 

workers recognized this feature of VO complexes but have encountered 
difficulty in the quantitative interpretation of the ESR line 
shapes because of the slow tumbling aspects inherent in viscous 

4 

nematic liquid crystals. Our approach is to analyse the ESR 
line shapes by the stochastic Liouville method as developed by 



2 


PBF for slow tumbling in an anisotropic media. 

Vanadyl acetylacetonate (VOAcAc) was chosen as our vanadyl 
spin probe because it has excellent ESR spectroscopic properties 
and because it has very nearly axially symmetric magnetic para- 
meters. In addition, our previous detailed work with VOAcAc 
in a viscous isotropic medium minimized experimental problems. 
The nematic mixture Phase V from EM Laboratories was the first 
nematogen studied in this work because it is nematic at low 
temperatures and has a wide nematic temperature range. A second 
nematic liquid crystal, butyl-p- (p-ethoxyphenoxy carbonyl) phenyl 
carbonate (BEPC) , which is a single component liquid crystal, 
was studied, and the results obtained were compared to those 
for the Phase V mixture. 
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Experimental 

Degassed samples of nematic liquid crystals doped with 
vanadyl acetylacetonate (VOAcAc) were prepared by standard 
vacuum line techniques. VOAcAc was purchased from Alfa Products 
and recrystallized from acetone. The nematic liquid crystal 
butyl-p- (p-ethoxyphenoxy carbonyl) phenyl carbonate (BEPC) was 
obtained from Eastman Kodak and recrystallized from methanol. 

The nematic mixture Phase V was purchased from EM Laboratories 
and used without further purification. VOAcAc concentrations 

_ I; 

(,5 x 10” M) were low enough so that contributions to the line- 
shapes from intermolecular interactions were negligible. 

ESR spectra were taken at X band with a Varian (V-4502) 
spectrometer system. 100 k Hz field modulation was used with low 
enough modulation amplitude that no distortion in lineshapes was 
observed. Microwave powers were low enough to preclude saturation 
effects. Temperature was controlled by a Varian E-257 variable 
temperature unit and measured with a copper-constantine thermo- 
couple. The temperature gradient over the active region of the 
cavity was 1-2°C. Spectral calibrations were obtained using a 
Brucker B-NM-12 gaussmeter or from the hyperfine splitting of 
VOAcAc in toluene at 2°C. 
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Theory 

PDF have developed a method for calculating slow tumbling 

ESR line shapes in an anisotropic solvent such as a nematic 
5 

liquid crystal. This method is based upon the stochastic 

Liouville method and is a generalization of the approach of 

7 

Nordio and co-workers for the motionally narrowed ESR line 

widths of a Brownian particle in liquid crystal. Here, we apply 

this general theory to vanadyl spin probes in nematic liquid 

crystals. Previous studies have given the appropriate spin 

3 8 

Hamiltonian for VOAcAc. ' In particular, the use of axially 

symmetric magnetic tensor parameters has been established. The 

_ 1 1 

motionally narrowed line widths (t < 8 x 10 sec/rad) calculated 
with axially symmetric magnetic parameters differ by less than 
0.5% from those calculated with completely asymmetric tensor 

O 

values. Alsu, our slow tumbling study in an isotropic solvent 

showed that the use oi axially symmetric: magnetic parameters gave 

excellent agreement with experimental line shapes for slow 

-11 -9 

tumbling t - 8 x 10 -1 x 10 sec/rad. Only with the appearance 

of features in the ESR line shapes that are characteristic of a 
rigid limit line shape does the use of completely asymmetric 
parameters improve the fit to experimental line shapes."^ Because 
VOAcAc aligns itself with the VO bond perpendicular to the 
director of the nematic liquid crystal, the rotation about the 
VO bond should be relatively fast and thus average the x and y axis 
magnetic asymmetry. 

The rotational diffusion for a molecule undergoing Brownian 
rotation in the presence of a cylindrically symmetric restoring 
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torque T is described by" 


H (n,t) - Rv2p < n 't> - ESIHT It CsinB(T)p(n,t>] 


where R is rotational diffusion coefficient, V 2 is the rotational 

diffusion operator, and 8 is the angle between the molecular Z~axis 

(VO bond) and the director. The director is taken to be parallel to 

the large dc magnetic field. For the axially symmetric VOAcAc, the 

g 

simple Maier-Saupe type potential 

U(B) = y cos 2 e (2) 

2 

is suitable and yields 

aft 

T = - ~ = 2 y sinScosS (3) 

op 2 

Note that for U(B)=0, Eq. (1) describes the familiar isotropic 

Brownian rotational diffusion. VOAcAc, an oblate symmetric top 

molecule, has 8 =tt/2 so that y is positive, 

eq 2 

From Eqs. (l)-(3), the equations for a nematic liquid crystalline 

5 

solvent may be obtained as in PBF by adding to the left hand side 
of each equation for the isotropic liquid [Eqs. (A.i)-(A.3) in 
reference 3]. 




G 


and by changing the right hand side of Eq. (A. 1) of reference 3 
to be 


< 

[ (2L+1)^ /I Q ] ^0^ o (O,&,O)exp[i x cos 3]sin8dB 


(5) 


where 


-r 

>■ 5 expt i ’■ 


cos B]sin8dB 


( 6 ) 


and X = -X 2 /kT. The absorption line shape Z is now given by 
8 ''if* 

Z" = J < 2L+1) ^ fif q 0 (S)exp[| x cos 8]sin8d8]cj^ Q ( j) 

j=lL even 0 ' 

where j refers to the eight allowed transitions. 

For anisotropic viscosity# 

2 L 


1 T, 




_ 1 _ 1 

in Eqs. (A. 1) - (A. 3) of reference 3 and x = xj_ in Eq. (4) above, 
x || is the rotational correlation time about an axis parallel to 
the director and xj_ refers to an axis perpendicular to the director, 


The equilibrium distribution Pq (£}) is given by 


22 r 2 

P n (fi) = exp(X cos 8)/4 tt ^exp(X cos 8)sinBd8 


( 8 ) 


The order parameter S is then 

S = <P 2 (B)> =^P 2 (8)P Q (n) dn 
where 

1 2 

P (8) “ 4(3 cos 8-1) 

2 1 


(9) 


(10) 
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Because VOAcAc is a much smaller molecule than the liquid 

crystal solvent molecules, the above diffusion equation for a 

Brownian particle may be expected to be modified for inertial 

effects. 10 One worthwhile modification to these equations involved 

the local solvent structure inducing a slowly flucuating torque 

11-13 

on our spin probe. 

The resulting equations explicitly contain inertial effects 

to lowest order and correspond to a model in which this torque 

flucuates between two values with a lifetime of t . These added 

m 

12 

equations are of the form 


[ (w-io, )- i m t ]C« „ _ _ + 
L u m J 0 , n , p ,m 


[L(L+1) [v/p+l cj f 1( „ ♦ 


l/p " X C 0 ,n,p-l,m] + [y- 


| 1 ] % 

T m 


C 


L 

0,n,p-l,m+l 


+ 


^ V™ C 0,n,p+X,m-1 ] = 


D . 0 


(ID 


where m=0, 1 and p=0, 1. I is the moment of inertia of the spin 
probe. Also 

T m = T T r. ( l/^T “D 

m L,n v L,n 


( 12 ) 
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Results and Discussion 

In order to simulate slow tumbling spectra, one must know 

the components of the g and A tensors which occur in the spin 

Hamiltonian. These values may vary slightly with solvent. The 

A 1 / and g- ft components are readily determined from a rigid limit 

spectrum. This spectrum was obtained by slowly cooling the liquid 

crystal sample through the nematic range with the dc magnetic 

field turned off and then taking the spectrum at 77 °k. The Ai 

and gx values may then be obtained from the isotropic ag and gg 

values gotten from the liquid crystal samples above the nematic- 

isotropic clearing point transition temperature. Our values for 

the magnetic parameters were essentially the same as those obtained 
3 8 

in toluene. ' For Phase V, Aj; = -184 G, Ax “ -68 G, g»i = 1.94 3, 
gx = 1.982. For BEPC, A // = -185 G, Ax® -68G, g u = 1.943, and 
gx = 1.982. 

Experimental ESR spectra were obtained for the nematic range 
72° to -29°C for Phase V and 82° to 37°C for BEPC. All spectra 
contained very asymmetric line shapes indicative of the slow 
tumbling region. In fact, the isotropic spectra taken just above 
the nematic-isotropic clearing point were obviously in the slow 
tumbling region. Although there was considerable supercooling 
below the nematic-solid freezing point, no discontinuities were 
observed in the ESR spectra or in their subsequent analysis. 

Simulated spectra corresponding to a Brownian particle in a 
liquid crystal solvent were calculated from Eq. (7) and Eqs. (A.l)- 
(A. 3) of reference 3 modified by Eq. (B.3) of reference 3 for non- 
secular effects and by Eqs. (4) and (5) for an anisotropic solvent. 
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Spectra were simulated by varying essentially only two parameters: 
t , the rotational correlation time, and A, which is proportional 
to the orienting potential. In general, the effect of increasing 
/A/ {A is negative for VOAcAc)is to decrease the separation between 
hyperfine lines and to decrease the line widths. An increase in 
t will generally broaden lines, decrease line separation (in the 
slow tumbling region) and increase asymmetry in the lines. Although 
the simulated spectra were in reasonable agreement, chere were 
significant and consistent descrepancies between calculated and 
experimental spectra. The calculated high-field "hyperfine lines" 
are more narrow than observed experimentally while the low-field 
lines were too broad. When Eq. (8) for anisotropic viscosity was 
incorporated into the above equations, the simulated spectra were 
in very good agreement with the experimental spectra. The spectra 
were in this case fit by varying three parameters x„ , t a , and A. 

It was noted that as the temperature was lowered the ratio 
of tj.Ah increased. Also, the best fits indicated that tji was 
becoming shorter with decreasing temperature. Because these results 
seemed implausible, the flucuating torque approach was tried. 

This analysis, although approximate, attempts to account for the 
small size of spin probe relative to the liquid crystal molecular 
size. Here we now incorporate Eqs. (12) and (13) into the above 
Brownian rotation equations. Very good fits were obtained by 
letting C L( „-2 and C L>n = 
t and A. See Table I. 


10(n^0) and varying two parameters. 
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TABLE la 

Ordering Parameter (S) and Rotational Correlation Time (t) 
as a Function of Temperature for Phase V 

Tomp (°C) S a 'k T a (xl0" scc/rad) 


69 

-0.11 

.11 

55 

-.14 

.20 

41 

-.18 

CO 

CO 

• 

34 

-.22 

.74 

20 

-.26 

1.4 

6 

-.29 

2.3 

-9 

-.33 

4.7 

-26 

-.36 

9.2 


a) Obtained from simulated jpectra using flucuating torque 
analysis with £- ft =2 and £ T =10 (nj^O) . 

Ju/ U § II 

b) Calculated numerically from Eq. (10) . 
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Table lb 


Ordering Parameter (S) and Rotational Correlation 
Time ( t ) as a Function of Temperature for BEPC 


Temp (°C) 

s a,b 

a -9 

T 3 ( Xl 0 

81 

-.16 

0.18 

75 

-.18 

0.25 

65 

-.20 

0.34 

56 

-.22 

0.57 

46 

-.25 

0.91 

38 

-.28 

1.2 


sec/rad) 


a. BEPC liquid crystal. Obtained from simulated spectra using 
flucuating torque analysis with £ T n *2 ar i =10 (n^O). 

u Lt § n 

b. Calculated numerically from Eq. (10) . 
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Conclusion 


The stochastic Liouville method as developed by PBF has been 
applied to the calculation of slowly tumbling ESR line shapes for 
vanadyl complexes in nematic liquid crystals. Quantitative 
information about the ordering and rotational diffusion was 
obtained for VOAcAc in Phase V for the nematic temperature range 
69°C (S = -0.11, t - 1.1 x 10 5|§) to -26°C (S- -.36, 


= 9.2 x 10 


_ 9 


sec 

rad 


rad 

) . For BEPC, the nematic temperature range 

-10 


sec 


rad 


81°C to 38°C corresponded to values (S = -0.16, x = 1.8 x 10 
and(s = -0.28, x = 1.2 x 10 |~j)respectively. Deviations from 
Brownian rotational diffusion was noticed and unquestionably arose 
from the relatively small size of the vanadyl probe. Application 
of a flucuating torque model approximated this effect and improved 
the fits to experimental spectra. 


The success of this approach encourages us to investigate 


larger VO probes (which would more nearly minic the liquid crystal 
itself) as well as to pursue the study of the highly viscous 
smectic liquid crystals. 

Also, the quantitative analysis of the spin probe ordering 
(and rotational motion) would allow one to monitor the ordering 
of liquid crystal subjected to both static and alternating 


electric fields. 
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